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ABSTRACT 

Spatially resolved polarized millimeter/submillimeter emission has been observed in the 
disk of HL Tau and two other young stellar objects. It is usually interpreted as coining from 
magnetically aligned grains, but can also be produced by dust scattering, as demonstrated ex¬ 
plicitly by Kataoka et al. for face-on disks. We extend their work by including the polarization 
induced by disk inclination with respect to the line of sight. Using a physically motivated, 
semi-analytic model, we show that the polarization fraction of the scattered light increases 
with the inclination angle /, reaching 1/3 for edge-on disks. The inclination-induced polar¬ 
ization can easily dominate that intrinsic to the disk in the face-on view. It provides a natural 
explanation for the two main features of the polarization pattern observed in the tilted disk 
of HL Tau (/ ~ 45°): the polarized intensity concentrating in a region elongated more or less 
along the major axis, and polarization in this region roughly parallel to the minor axis. This 
broad agreement provides support to dust scattering as a viable mechanism for producing, at 
least in part, polarized millimeter radiation. In order to produce polarization at the observed 
level (~ 1%), the scattering grains must have grown to a maximum size of tens of microns. 
However, such grains may be too small to produce the opacity spectral index of < \ ob¬ 
served in HL Tau and other sources; another population of larger, millimeter/centimeter-sized, 
grains may be needed to explain the bulk of the unpolarized continuum emission. 
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1 INTRODUCTION 

Polarized millimeter/sub-millimeter emission has been observed in 
the disks around 3 young stellar objects: IRAS 16293-2422B (Rao 
et al. 2014), HL Tau (Stephens et al. 2014) and L1527 (Segura- 
Cox et al. 2015). It is usually interpreted as coming from magnet¬ 
ically aligned dust grains (e.g., Cho & Lazarian 2007). This inter¬ 
pretation appears consistent with the data on IRAS 16293-2422B 
and LI527, where the observed polarization patterns are broadly 
similar to those expected from grains aligned by a predominantly 
toroidal magnetic field (Rao et al. 2014; Segura-Cox et al. 2015; 
see § 5 for more discussion). In contrast, the polarization vectors 
in the disk of HL Tau (which is inclined with respect to the line 
of sight by ~ 45°, ALMA Partnership: Brogan et al. 2015; Kwon 
et al. 2011) are all roughly parallel to the minor axis (see Fig.l), 
which is not compatible with the pattern expected of grains aligned 
by a toroidal disk field, although more complicated field configu¬ 
rations cannot be ruled out (Stephens et al. 2014). Another draw¬ 
back of the toroidal field-aligned grain model is that it predicts a 
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lower polarization fraction along the major axis than along the mi¬ 
nor axis (see Cho & Lazarian 2007, Fig. 10), which is the opposite 
of what is observed (Fig. 1). A much better fit is provided by a uni¬ 
directional magnetic field approximately (within ~ 10°) along the 
disk major axis. However, such a field configuration would be diffi¬ 
cult to maintain against the disk differential rotation. Furthermore, 
there is growing evidence for grain growth in protoplanetary disks, 
up to millimeter or even centimeter sizes (e.g., Perez et al. 2012; 
Alecian 2013; Testi et al. 2014). It is unclear whether such large 
grains would be aligned with respect to the magnetic field through 
the currently favored mechanism of radiative torque because their 
magnetic moments may not be large enough to provide the fast pre¬ 
cession needed (Lazarian 2007) and their slow internal relaxation 
makes the alignment less efficient (Hoang & Lazarian 2009). These 
difficulties provide a motivation to investigate other potential mech¬ 
anisms for producing, at least in part, the polarized mm emission 
from the HL Tau disk. 

One possibility, first investigated in detail by Kataoka et al. 
(2015a), is that the polarized millimeter disk emission comes from 
scattering of large dust grains. Kataoka et al. stressed the need 
for anisotropic radiation field for this mechanism to work effi- 
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Figure 1. Observed 1.3 mm polarization pattern of the HL Tan disk. Plotted 
are total intensity (contours), polarized intensity (color map, > 2cr), and 
polarization vectors (> 3cr), with length proportional to the polarization 
fraction. Adapted from Stephens et al. (2014). The major axis determined 
from the ALMA data (ALMA Partnership: Brogan et al. 2015) is tilted by 
10° clockwise with respect to the vertical direction. 


ciently. As examples of this requirement, they considered radiation 
from structured disks with either an axisymmetric ring or a non- 
axisymmetric lobe, both viewed face-on. Here, we extend Kataoka 
et al.’s work to show that significant polarization is produced even 
in a smooth disk, as long as the disk axis is tilted away from the line 
of sight by a large enough angle. The physical reason can be under¬ 
stood most easily in the limits of optically and geometrically thin 
dust emission and Rayleigh scattering, where the degree of polar¬ 
ization is peaked near a scattering angle of 90°. In these limits, the 
radiation field to be scattered by dust grains becomes essentially 
two dimensional, concentrated in the disk plane. In such a case, 
only the radiation propagating along the disk major axis would be 
scattered by 90° to reach the observer and be maximally polarized, 
with a polarization direction along the minor axis in the plane of 
the sky. This simple geometric effect provides a natural explana¬ 
tion for the two main features of the polarization pattern observed 
in the HL Tau disk (see Fig. 1): (1) polarized emission concentrat¬ 
ing in an elongated region more or less along the major axis, and 
(2) polarization vectors in this region roughly parallel to the mi¬ 
nor axis. These are the features that are difficult to explain in the 
current models of polarization from magnetically aligned grains ^ 
They provide a strong motivation to explore in more detail the po¬ 
larization pattern produced by dust scattering. It is a first step to¬ 
ward disentangling the contributions of dust scattering and grain 
alignment to the disk polarization, which is needed in order to take 
full advantage of the high resolution ALMA and JVLA polariza¬ 
tion observations that will become available soon for probing the 
grain growth and/or magnetic field in protoplanetary disks. 

In the remainder of the paper, we will put the above qual- 


^ Grains aligned by a toroidal magnetic field can explain the orientations 
of the polarization vectors near the center, but not at larger distances (see 
the left panel of Fig. 2 of Stephens et al. 2014). 


itative picture on a more quantitative ground by computing the 
scattered radiation from an inclined disk (such as the HL Tau 
disk) semi-analytically using the so-called “thin-disk” approxima¬ 
tion that brings out the essential physics transparently (§ 2); a more 
complete treatment of the radiative transfer problem will be post¬ 
poned to a future investigation. The effects of disk inclination on 
the polarization pattern will be emphasized. The scattered radia¬ 
tion is compared with the direct dust emission in § 3, where the 
fractional polarization of the total intensity is computed. In § 4, we 
comment on the implications of the polarization fraction observed 
in HL Tau on the grain size distribution. We discuss the possibility 
of dust scattering-induced polarization for other sources and future 
model refinements in § 5 and conclude in § 6. 


2 POLARIZED RADIATION FROM DUST SCATTERING 

In this section, we develop a semi-analytic model for polarization 
from dust scattering under the simplifying assumptions that the 
dust emission is optically thin and that the dust disk is geometri¬ 
cally thin, with a local scale height much less than the radius. The 
latter would be a particularly good approximation if the relatively 
large grains responsible for the emission and scattering of millime¬ 
ter radiation are settled toward the disk mid-plane. To isolate the 
polarization due to dust scattering from that caused by grain align¬ 
ment, we further assume that the grains are not aligned, so that the 
radiation directly emitted by dust is unpolarized. Under these as¬ 
sumptions, the incoming radiation field seen by the scattering dust 
grains at a given point inside the disk can be decomposed into two 
components: a local, roughly isotropic, component emitted by the 
(non-aligned) grains in a region of size comparable to the (dust) 
scale height surrounding the point, and an anisotropic component 
coming from the grains in the rest of the disk beyond the local 
region. The former does not lead to significant polarization after 
scattering because scattering of isotropic radiation by non-aligned 
grains does not have any preferred direction. It is the latter that is 
mainly responsible for the polarized radiation after scattering. The 
computation of this component will be drastically simplified by the 
thin-disk approximation, as we show below. 


2.1 Formulation of the Problem 

To facilitate the computation of the scattered light from dust grains 
at a location r inside a disk with axis inclined by an angle i with 
respect to the line of sight, we define two coordinate systems, both 
centered on the scatterer’s location r. The first is fixed in the ob¬ 
server’s frame, with horizontal x'- and vertical y'- axes in the plane 
of the sky, and z'-axis pointing toward the observer. The second is 
fixed on the disk, with the y-axis coinciding with the y'-axis, and the 
z- and x-axis rotating around the y'-axis counter-clockwise (viewed 
along the minus y'-axis) from the z'- and v'-axis, respectively, by 
the inclination angle i. We align the major axis of the disk in the 
plane of the sky along the y- (and y'-) axis (see Fig. 2 below for an 
illustration), and the disk normal direction along the z-axis. 

Since the scattered radiation is linearly polarized, we will 
compute the three Stokes parameters 7^, Q and U defined in the 
observer’s frame (v',y',z') separately, starting with the intensity Is 
(we drop the subscript v for Q and U since, in our model, they are 
assumed to be produced only by the scattered radiation, without 
any contribution from the direct radiation, unlike the total intensity 
/). In the optically thin limit, the intensity of the scattered radi¬ 
ation along the line of sight Is is given by an integration of the 
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source function of the scattered radiation, S s,z'y along the z'-axis. 
The source function S s,z' at any location r inside the disk is given 
by 

SsAr)= — f ^nr,e,4>)dil (1) 

o-s J aQ 

where /(r, 0, 0) is the intensity of the unpolarized incoming radi¬ 
ation seen by the scattering dust grains at the location r along the 
direction of polar angle 0 (measured from the z-axis or the disk nor¬ 
mal direction) and azimuthal angle 0 (measured counter-clockwise 
from the v-axis on the disk), cr^ is the solid angle Q-integrated (to¬ 
tal) scattering cross section, and dcrldCl is the differential cross sec¬ 
tion for scattering the radiation /(r, 6, 0) into the line of sight (i.e., 
along the z'-axis). 

The intensity /(r, 6, 0) is given by an integration of the source 
function for thermal dust emission, the Planck function By(T), over 
the optical depth dzats = KabsP dl (where Kabs is the thermal dust 
absorption opacity and p the mass density at a source location rO 
along a line in the direction {6, 0) up to the scatterer at r (the quan¬ 
tity / is the distance between ri and r): 

= Kabs(r\y^^-!^^^p{r\)dl, (2) 

where we have assumed that the photon energy hv is substantially 
less than kT so that the Rayleigh-Jeans law By{T) - 2v^kT /c^ is 
applicable; we have checked that, for the HL Tau model to be dis¬ 
cussed below in § 2.3, the results will not change significantly if 
the Planck function is used instead. 

Substituting the above expression for /(r, 6, 0) into Eq.(l) and 
reorder the integrals, we can rewrite the source function for the 
scattered radiation into: 

2v^k r“ r da- 

Ss.z' =-2 — I #1 dl \ de—Kabs(r\)p(r\)T(ri)sme 

c^o-, Jo Jo Jo da ( 3 ) 

=S 0 + S oa, 
where 

2y^k r do- 

Sq =-— <i0 dl \ de—Kabs(Xx)p{Yx)T{Yx)^me, (4) 

cV, Jo Jo Jo dkl 

and 

2v^k C r do- 

Soo=—. J0 dl \ de — Kabs(yx)p{Yx)T{Yx)sme. (5) 

cV, Jo Jh Jo dkl 

The two quantities, Sq and S oo, denote the contributions to the 
unpolarized incoming radiation to be scattered at the location r 
from two conceptually distinct regions respectively. For a geometri¬ 
cally thin dust disk, the unpolarized radiation coming from a region 
within a distance on the order of the local (dust) scale height H is 
expected to be more or less isotropic. This near-field contribution, 
denoted by ^o, produces little polarized radiation after scattering; it 
will be discussed in the next section, together with the unpolarized 
direct dust emission. In contrast, the thermal dust emission coming 
from well outside this local region (i.e., I » H) is mostly con¬ 
fined close to the disk plane. This far-field contribution, denoted 
by 5oo, is highly anisotropic, with the bulk of the radiation beamed 
into a narrow range of polar angle near 6-7:12 along any az¬ 
imuthal direction 0. Specifically, the vertical column of disk mate¬ 
rial passing through the source location ri at a distance / from the 
scatterer contributes radiation only within a range of polar angle 
66 ~ 2H{yx)II « 1, where H{yx) is the (dust) scale height at ri. Re¬ 
placing the integral over angle 6 in the expression for Soo (eq. (5)) 
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by this rough estimate and making use of6 ^ 7:12, we have approx¬ 
imately 


2y^kKabs 


f 


, do- . , 

(p), 

dkl 


( 6 ) 


where we have assumed a spatially constant absorption opacity Kabs 
for simplicity. The auxiliary quantity A(r, 0) is a line integral along 
the direction of constant 0 in the disk plane defined as 


A(r, 0) 


^(X 

Jh 


dl 


^r i)r(ri) 
/ 


(7) 


where E(ri) = 2p(yx)H(yx) is the column density at the source loca¬ 
tion ri. 

One can determine the distance between the source location ri 
and the center of the disk through 


ri - sjr^ + B - 2rl cos(0 - 0r), 


( 8 ) 


where r and 0r are the radius and azimuthal angle of the scatterer 
in an (x, y, z) coordinate system that centers on the star (rather than 
the scatterer at r). For the axisymmetric disk that we will consider 
below, the radius r/ uniquely determines the column density E and 
temperature T that appear in eq.(7). 

Once the integral A(r, 0) is computed, the only term that is left 
to determine in eq.(6) is the differential cross section for scattering 
Jcr/JO. For illustrative purposes, we will consider the dust scat¬ 
tering under the Rayleigh approximation, which is valid when the 
grain sizes are smaller than the wavelength divided by 27: (see, e.g.. 
Fig. 10 of Canovas et al. 2013 and Fig. 6 below); we will check in 
§ 4.1 that this condition is satisfied for the HL Tau model to be dis¬ 
cussed in the next section. In this limit, the differential cross section 
is given by 


Jcr 


^(i 

16;r ^ 


-h cos' 


00 


(9) 


where 6s is the scattering angle between the incoming radiation 
along the direction {6, 0) and the outgoing radiation along the line 
of sight (i.e., the z'-axis), given by 


cos 6s - sin i cos 0. 


( 10 ) 


After scattering, a fraction of the initially unpolarized incoming 
radiation becomes polarized. The polarization fraction is 


1 - cos^ 6s 
1 -h cos^ 6s ’ 


( 11 ) 


which peaks at - t:/2. 

The polarization direction is perpendicular to the scattering 
plane formed by the incoming direction (6, 0) and the z'-axis. It lies 
in the x'-y' plane of the sky, at an angle 0' -h 7:12 measured counter¬ 
clockwise from the v'-axis, with the angle 0' given by 


cos i cos 0 


^ sin^ 0 -h cos^ i cos^ 0 


( 12 ) 


This linearly polarized radiation is the source of the observed 
polarized radiation in the plane of the sky, through the source func¬ 
tions for the Stokes parameter Q and U : 

Iv^kKabs r^^j^daiej 

Sq,o^~ - 5 - d(p——A(r,(p)p(ejcos(2(p), (13) 

cVj Jo da 

and 

2v^kKabs C^jj.daiej . 

Sv.«>~ - 5 - d(p——K(r,<p)p{ejsm(2(p). (14) 

Jo da 
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For the optically thin radiation that we are considering, Q and U 
along the line of sight passing through any location r in the disk 
is simply given by their respective source functions, multiplied by 
the optical depth for scattering Ar^ Ksca^ix)! cos i (where Ksca is 
the scattering opacity) through the disk: Q - S q^oo^Ts and U - 
S f/^ooAr^. 


2.2 Inclination-Induced Polarization 


As a simple illustration of polarization from scattered light induced 
by disk inclination, we consider the limiting case where the incom¬ 
ing radiation field seen by the scatterer at the location r is confined 
to an infinitely thin disk plane (so that the near-field contribution 
to the scattering source function, Sq, can be ignored compared to 
the far-field contribution Soo) and is isotropic in the azimuthal (0) 
direction (but highly anisotropic in polar angle 6). In this case, the 
integral A(r, 0) defined in eq. (7) is independent of 0 and can be 
moved outside the integral over 0 in the source functions for the 
total scattered intensity (eq. (6)), Q (eq. (13)) and U (eq. (14)), 
so that 


, 2v^kK,l,,A f"" ,, 3 . 2 . 2 

Soo- -- d(p^—{l + sm icos 0) 

Jo 167r 

J0(1 -h sin^ i cos^ 0) = 7tC(2 + sin^ i), 

0 


(15) 


where C = Sv^kKabs-^/i^TTc^) is a constant independent of the an¬ 
gles i and 0, and 


> 2,00 


-c 


/ 


J0(1 - sin /cos 0)- 


'C0s^/-hsin 0 (16) 


: 7tC sim /, 


> u,oo 


-c 


f 


J0(l-sin^ /cos^ 0) 


cos / sin(20) 
cos^ 0 cos^ / -h sin^ 0 


: 0. (17) 


The fact that S u,oo is zero and S q^oo is positive (for / ^ 0) means 
that the inclination-induced polarization is always along the v'-axis 
(or the minor axis of the disk) in the plane of the sky for Rayleigh 
scattering (see Fig. 2 for the x'-y' coordinates, and eq. (23) below 
for the relation between the Stokes parameters Q and C/, and po¬ 
larization angle a). This is expected physically because, in a tilted 
disk, the light coming from different directions in the disk plane 
will be scattered by different angles toward the observer. In partic¬ 
ular, the light coming from a direction along the major axis will 
always be scattered by njl. In the Rayleigh limit, this part of the 
light will be fully polarized along the minor axis of the disk. In con¬ 
trast, the light coming from a direction along the minor axis will be 
scattered by either njl - i or n/l + i. This part of the light will be 
partially polarized along the major axis, with a polarization fraction 
of cos^ //(I -h sin^ /). The difference in the fraction of polarization 
leads to more scattered light polarized along the minor axis than 
along the major axis, despite the fact that the scattering cross sec¬ 
tion is less for the former than the latter (see eq. (9)). This generic 
tendency for the inclination-induced polarization to align with the 
minor axis provides a natural explanation for the polarization di¬ 
rections observed in HL Tau (see Fig. 1). 

It is easy to determine the fraction of the total scattered light 
that is polarized along the minor axis: 


Ps(i) 


h 


S 


2 


- 

^ U,oo 


Soo 


sin^ / 

2 -h sin^ / 


(18) 


This same expression can be obtained if one considers only the ra¬ 
diation coming from directions along the major and minor axes. 
Note that the maximum degree of polarization reaches 1/3 when 
the disk is viewed edge-on (/ = 90°). For the inclination angle 
i - 0°, 30°, 45° and 60° to be considered in the next subsection, 
the fractional polarization is Ps = 0, 1/9, 1/5, and 3/11, respec¬ 
tively. The analytically obtained polarization fraction in this simple 
limiting case will be used to interpret the results obtained numeri¬ 
cally in more general cases. 


2.3 Intrinsic Polarization from Azimuthally Anisotropic 
Radiation: an Example 

As emphasized by Kataoka et al. (2015a), the radiation field in the 
disk plane is not isotropic in general, and the anisotropy in the az¬ 
imuthal (0) direction leads to polarized scattered light even in the 
face-on case. In a tilted disk, the observed polarization pattern is 
expected to be shaped by the interplay between those produced 
by anisotropy in 0-direction and disk inclination (which relies on 
strong anisotropy in 6 direction). To illustrate this interplay, we 
adopt Kwon et al. (2011) model of the HL Tau disk, where the 
distributions of temperature and column density as a function of 
the cylindrical radius R are parametrized as 

(19) 


— 


3/2-p-q/2 


exp 


Il2-p-ql2 


( 20 ) 


where Rc is a characteristic disk radius beyond which the column 
density drops off exponentially, and p is an exponent that, together 
with the exponent q, controls the column density distribution; it is 
not to be confused with the polarization fraction. The temperature 
profile yields a thermal scale height for the gas 


H^Ho 



j3/2-^/2 


( 21 ) 


Although higher resolution ALMA observations have revealed sub¬ 
structures (rings) on the disk (ALMA Partnership: Brogan et al. 
2015), the above is still the best model at the CARMA resolu¬ 
tion that was used to detect the polarization in HL Tau. Adopting 
q - 0.43, Kwon et al. found a set of parameters that best fit their 
CARMA observations at 1.3 and 2.7 mm: Rc ^ 79 AU, ~ 1, and 
a scale height for gas at Rc of Ho ~ 16.8 AU (see also Kwon et al. 
2015). Since the relatively large grains responsible for the scattered 
radiation may settle toward the midplane, the scaling factor Ho in 
eq.(21) may need to be reduced by some (potentially large) factor 
(Kwon et al. 2011). We have experimented with reduction factors 
of 1, 10, 50 and 100, and found very similar results. In what fol¬ 
lows, we will focus on the case where the scale height is the same 
for the dust and gas. 

With the disk structure specified, we can now compute the 
Stokes parameters Q, U and from their source functions given 
by eqs. (13), (14) and (6). In this section, we consider the contribu¬ 
tion A oo to the total observed intensity from the far-field scattering 
source function S oo only, in order to facilitate comparison with the 
analytic results obtained in the preceding subsection; the contribu¬ 
tions from S 0 and direct dust emission would lower the polariza¬ 
tion fraction, and will be considered in the next section. From these 
three Stokes parameters, we can determine the total intensity of the 
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polarized radiation 

Ip = ( 22 ) 

the polarization angle a (measured counter-clockwise from the x'- 
axis in the plane of the sky) 

or = ^atan2|^|, (23) 


polarization and inclination-induced polarization. It provides a nat¬ 
ural explanation for the distribution of polarized intensity observed 
in HL Tau^, which has an inclination angle of i ~ 45° (ALMA 
Partnership: Brogan et al. 2015; Kwon et al. 2011). 

2.4 Interplay Between Intrinsic and Inclination-Induced 
Polarization 


where the function atan2 returns the appropriate quadrant of the 
computed angle based on the signs of Q and C/, and the polarization 
fraction 


Ps,oo = 



(24) 


To understand the interplay between the intrinsic polarization and 
inclination-induced polarization more quantitatively, we plot in 
Figure 3 the distribution of a dimensionless quantity 

PQ^-p- ( 25 ) 

^ 5,00 


which can be compared directly with that given analytically in 
eq. (18). 

In Fig.2, we show the polarization vectors and the spatial dis¬ 
tribution of the polarized intensity Ip for / = 0°, 30°, 45°, and 60°. 
The intensity is measured in units of TZlKabsKscay^kToI, where 
So and Tq are the characteristic column density and temperature 
of the disk, and Kabs and Ksca the absorption and scattering opac¬ 
ity. In the face-on (/ = 0°) case, the fraction of polarization of the 
scattered light is zero at the center because the light to be scat¬ 
tered there comes isotropically along all azimuthal directions for 
the prescribed axisymmetric disk. The radiation seen by the scat- 
terer becomes more and more beamed in the radial direction as the 
radius increases, because of the drop in temperature and column 
density. As a result, the light is polarized in the azimuthal direction 
and the polarization degree increases outward. We will refer to the 
polarization induced by anisotropic radiation in the azimuthal di¬ 
rection in the face-on case as the “intrinsic” polarization. Note that 
although the polarization fraction ps^oo can reach a value as high as 
50% or more near the outer edge, the total scattered polarized in¬ 
tensity Ip is rather low. As a result, the intrinsic polarization from 
dust scattering can be easily modified, indeed dominated, by the 
inclination-induced polarization. 

When the disk is tilted away from the line of sight, both the po¬ 
larized intensity and orientations of the polarization vectors change 
drastically compared to the face-on case. As Fig. 2 shows, the po¬ 
larized intensity peaks at a ring in the face-on case, with the in¬ 
ward decrease caused by radiation becoming more isotropic in the 
azimuthal direction and the outward decrease from the exponen¬ 
tial drop-off in column density. In contrast, the polarized inten¬ 
sity in the i = 30° case peaks in the central region, as a result 
of the inclination-induced polarization. The polarization vectors 
in this (central) region (within ~ 0.5Rc of the origin) lie more or 
less along the minor axis, consistent with the analytic results for 
the inclination-induced polarization derived in the last subsection. 
Outside the central region, the polarization directions are broadly 
similar to those in the face-on case, indicating that the intrinsic 
polarization remains important there. A difference is that the ax¬ 
isymmetric polarization pattern in the face-on case becomes highly 
non-axisymmetric in this moderately tilted case, with both the po¬ 
larized intensity and the polarization fraction significantly higher 
along the major axis than along the minor axis. 

As the disk tilt angle increases further, the polarization pat¬ 
tern becomes more dominated by that induced by inclination. Go¬ 
ing from i = 30° to 45° to 60°, we see a clear trend for increasing 
polarized intensity in the central region, a larger fraction of the po¬ 
larization vectors parallel to the minor axis, and more elongation of 
the polarized intensity along the major axis (see Fig.2). The elon¬ 
gation is a generic feature of the interplay between the intrinsic 


along the major and minor axes for the i = 45° case. Since U = 
0 along the major and minor axes, this quantity is essentially the 
polarization fraction ps,oo defined in eq.(24), except that it retains 
the sign of the Stokes parameter Q. A positive (negative) pq means 
that the polarization is along the v'-axis (y'-axis) in the plane of the 
sky. 

In the face-on case (/ = 0°), the intrinsic polarization frac¬ 
tion increases monotonically from zero to a large value approaching 
unity as the distance from the origin increases along the major (or 
y'-) axis, with the polarization vectors parallel to the v'-axis in the 
plane of the sky. The quantity Pq is thus positive, as shown in Fig. 3 
(solid black line). Along the minor axis, the fractional polarization 
is the same, but the polarization direction is along the y'-axis, with 
a negative Pq (dashed black line). For comparison, we plot in the 
same figure the numerically computed distribution of pq along the 
major and minor axes as a function of the (de-projected) distance 
from the origin for the i - 45° case. Clearly, the inclination has 
shifted the curves for the intrinsic (i = 0°) Pq along both major and 
minor axes upward, which is not surprising since inclination tends 
to produce polarization along the v'-direction (corresponding to a 
positive 2, see eq.l6). The amount of upward shift is different at 
different locations, however. In the central region where the intrin¬ 
sic polarization fraction is relatively low (within about 0.5 R^), the 
shift is by ~ 0.2 along both the major and minor axes; it is the value 
expected from the analytic results for an azimuthally isotropic ra¬ 
diation field given in eq. (18). The behavior in the central region 
where the intrinsic polarization is weak is therefore easy to under¬ 
stand; it is dominated by the inclination-induced polarization. 

The behavior outside the central region where the intrinsic po¬ 
larization fraction is higher is more complicated. It can be repro¬ 
duced exactly, however, by the formula 

, -Pe.'" 

Pq — - 0 -'- 9 - 9 -’ 

2 + sin^ i (1 + sin^ //2)[1 + (1 - pqj„) sin^ i/l] 

along both the major and minor axes. Note that the first term on 
the right hand side is simply the inclination-induced polarization 
given by eq. (18), and the second term is the intrinsic polarization 
fraction pqj^ (the face-on case) modified by the inclination. This 
formula can be derived heuristically under the assumption that the 
radiation seen by the scatterer comes from only two directions: the 
X- and y-axes in the disk plane. 

^ The observed direction of elongation does not lie exactly along the major 
axis, but is offset by a small angle of ~ 10° (see Stephens et al. 2014 and 
Fig. 1). This offset is not explained in our model under thin-disk approxi¬ 
mation, and may require full 3D models and/or additional physics, such as 
grain alignment. 
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Figure 2. Effects of disk inclination on the polarization of scattered millimeter radiation. Shown are the intensity of the polarized radiation (color map, in 
units of 2LQKabs/<scay^kTolc^) and polarization vectors (line segments, with length proportional to the polarization fraction ps,oo defined in eq.(24)), for / = 0°, 
30°, 45° and 60°. As the disk becomes more tilted, both the polarized intensity and polarization directions become more dominated by those induced by 
inclination. The elongation of the polarized intensity along the major axis and orientations of the polarization vectors along the minor axis in the high intensity 
region in the inclined cases are broadly consistent with the observed pattern in HL Tau (see Fig. 1). The numbers 0.5, 1.0, 1.5 and 2.0 in each panel measure 
the de-projected distances from the center in units of the characteristic radius Rc. The x'-axis and y'-axis defined in the plane of the sky are shown in the 
lower-right panel for reference. 


In the limit where the intrinsic polarization fraction IpqjJ 
0, we have 


Pq 


sin^ i 
2 -f sin^ i 


(27) 


In the opposite limit where the intrinsic polarization fraction 
\PQ,in\ ^ 1, we need to consider the major and minor axes sepa¬ 
rately. Along the major axis where pq is positive, we have pQ^n 
1, which yields 


Pq 


(28) 


Along the minor axis where pq is negative, we have pqj^ -1, so 


that 
Pq - 


1 - sin^ i 
1 -f sin^ i ’ 


(29) 


polarization patterns along the major and minor axes. It becomes 
more extreme as the inclination angle i approaches 90°. 

Another difference is that there exists a point of zero polar¬ 
ization on the minor axis across which pq changes sign (or the 
polarization direction changes by 90°). This null point occurs at a 
location where 


PQJfi 


2 - sin i 


(30) 


which has a value of -1/3 for i - 45°. Indeed, if the inclination 
angle i is known independently, one can in principle deduce the 
intrinsic value of Pq along the major and minor axes from the value 
of Pq through 


PQ,in - 


Pq - sin^ //(2 -f sin^ i) 

[(2 - sin^ i) + Pq sin^ /]/(2 -f sin^ i) 


(31) 


which, for i - 45° shown in Fig. 3, approaches -1/3. The differ¬ 
ence in the asymptotic behavior of Pq, particularly the polarization 
fraction \pqI highlights one of the major differences between the 


However, it is difficult to infer the value of pq from observation di¬ 
rectly because it is the observed Stokes Q parameter normalized by 
the scattered intensity from the far-field, /^ oo, which cannot be mea- 
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Figure 3. Dimensionless quantity pQ defined in eq. (25). Plotted are the 
distributions of pQ along the major and minor axes, denoted by the sub¬ 
script “maj” and “min” respectively, as a function of de-projected distance 
from the origin for the face-on i = 0° case (solid and dashed black) and 
the i = 45° case (solid and dashed blue). Its polarization fraction \pq\ ap¬ 
proaches 20% in the central region and, at larger distances from the origin, 
is substantially larger along the major axis than along the minor axis. 


sured directly. What can be measured is the total intensity, which 
we discuss next. 


3 TOTAL INTENSITY AND POLARIZATION FRACTION 


The total intensity I of the radiation along the line of sight is the 
sum of the direct dust emission /j, the scattered radiation from near¬ 
field Isfi, and the scattered radiation from far-field Is,oo- The source 
function for the near-field contribution at a location r inside the 
disk can be estimated approximately assuming isotropic incoming 
radiation from within a uniform sphere of radius H, the local scale 
height: 


v^te„i,jS(r)7’(r) 


>0 ■ 


(32) 


which is multiplied by the scattering optical depth Ar^ 
Ksca^ir)/ COS i to yield the intensity 


I 


V^kKabsKsca^^(r)T{Y) 


s,0 


(33) 


For the HL Tau disk model discussed in § 2.3, the intensity of the 
scattered radiation A o from the near-field is weaker than that from 
the far-field oo everywhere for i - 45°, as can be seen from Fig. 4, 
where the ratio of the two is plotted. The ratio peaks in a ring be¬ 
tween ~ 0.5/?c and ~ IRc (where Rc is the characteristic radius of 
the disk), with a maximum value of ~ 60%. Near the peak, the 
unpolarized scattering intensity Is,o reduces the polarization frac¬ 
tion by a factor of ~ 1.6, from ~ 20% to ~ 12%. Going outward 
from the ring, the ratio drops rapidly because the near-field inten¬ 
sity Is,o is determined by the local column density, which drops off 
exponentially with radius, whereas the far-field intensity Is,oo is de¬ 
termined globally, including contributions from the bright central 
region that decrease with radius more slowly than exponential. In¬ 
side the ring, the ratio decreases with decreasing radius because of 
a smaller scale height H, which decreases the size of the region 


Figure 4. Ratio of the near-field and far-field contributions to the intensity 
of the scattered radiation, A,o/A,oo, for the i = 45° case. The far-field contri¬ 
bution Is,oo dominates the near-field contribution Isfi everywhere, especially 
in the central region and near the outer edge. 


where the incoming radiation for the near-field source function S o 
comes from relative to that for the far-field source function Aoo. In 
any case, the polarization fraction of the scattered radiation remains 
high, of order 10% or more, after both the near- and far-field con¬ 
tributions are taken into account. It is much higher than observed 
in HL Tau (of order 1%), and needs to be further reduced, by the 
unpolarized direct dust emission^. 

The source function for the direct emission is simply the 
Planck function By(T) ^ Iv^kTjc^. The optical depth for dust ab¬ 
sorption along the line of sight passing through a location r in the 
disk is At« ^ Kabs'^i^)! cos i under the thin disk approximation. To¬ 
gether, they yield the intensity for the direct emission 

^ 2y^kKabs^(v)T(v) 

fd ~ ^^-• (34) 

cos I 

This estimate allows us to evaluate the ratio of the intensities of the 
scattered and direct emission: 


h 

Id 


ls,0 "f Is, 


1 3 r^’ 

2 ^ 167r Jo 


V 

' ^sca^O ^ ^ 

2-0 

r 

d(p{l + sin^ /cos^ 0) 


p dl 

Jb I 


^rO r(r0 


E(r) T{v) 


(35) 


Note that the ratio E(r)/Eo and the second term inside the square 
bracket are dimensionless quantities that depend only on the shape 
of the column density and temperature profiles. The overall scaling 
is set by the characteristic scattering optical depth Ar^^c = Ksca^o- 
In order to reduce the high polarization fraction of the scattered 
radiation at 1.3 mm to the observed value of about 1%, we need a 
rather small value of Ar^^c ~ 0.07, so that the scattered radiation is 
heavily diluted by the unpolarized direct emission. 

We stress that the inclusion of the unpolarized radiation A o 
and Id changes neither the polarized intensity nor the polarization 
direction shown in Fig. 2 (the lower-left panel). What is changed is 
the polarization fraction. In Fig. 5, we plot the distribution of the 


^ Direct emission can be polarized, due to, for example, grain alignment. 
Polarized direct emission needs to be included in a more complete model 
(as discussed in § 5). 
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Figure 5. Model for the polarization of HL Tau disk. Plotted are the polar¬ 
ized intensity (color map, in units of 2L^KabsKsca^^kTQlc^) and polarization 
vectors (line segments, with length proportional to the polarization fraction 
of the total intensity) for i = 45°. A characteristic scattering optical depth 
Ar^^c ~ 0.07 is needed to bring the high polarization fraction of tens of 
percent relative to the scattered intensity (shown in the lower-left panel of 
Fig. 2) down to the observed level of ~ 1% relative to the total intensity. 

polarized intensity for the i - 45° case, as in the lower-left panel 
of Fig. 2, but with the length of the overlaid polarization vectors 
scaled by the new polarization fraction (relative to the total inten¬ 
sity). This figure represents our final model for the FIL Tau disk. It 
has three features that are broadly consistent with observations: (1) 
the region of high polarized intensity is elongated along the major 
axis, (2) the polarization vectors in this region are nearly parallel 
to the minor axis, and (3) the polarization fraction in the region is 
about 1%. Along the ridge of detectable polarized intensity, the ob¬ 
served polarization fraction appears to be somewhat higher toward 
the edge of the disk (at a distance of ~ ~ 80 AU) than near the 

center, although it is unclear how significant the trend is in view of 
the low polarized intensity near the edge. This trend was present in 
our model when only the scattered radiation was considered (see 
Fig. 2 and 3), but was washed out by the total intensity in Fig. 5. 

The semi-analytic theory that we have developed so far under 
the assumption of geometrically thin disk, optically thin emission, 
and only Rayleigh scattering is independent of the detailed proper¬ 
ties of dust grains. This independence makes the broad agreement 
between our model and the main polarization features observed in 
HL Tau rather robust. In the next section, we will try to put con¬ 
straints on the grain size distribution in the HL Tau disk, which is 
much more uncertain. 


4 IMPLICATIONS ON GRAIN SIZE 

4.1 Scattering Opacity and the Need for Large Grains 

The main free parameter of the final model for the HL Tau disk 
polarization discussed in § 3 is the characteristic scattering op¬ 
tical depth which controls the polarization fraction. In or¬ 

der to produce the observed polarization fraction of ~ 1%, a 
value of ^ 0.07 is required. This value yields a scattering 
opacity Ksca ~ 10“^ cm^ g“^ (cross section per unit total, rather 
than dust, mass) at 1.3 mm (the wavelength of the HL Tau disk 


polarization observation), using the characteristic column density 
Eo = 68 g cm“^ from the best-fit disk mass of 0.13 M© of Kwon et 
al. (2011). This scattering opacity can put constraints on the grain 
size distribution, although they depend on the dust composition, 
which is uncertain. As an illustration, we consider the model of dust 
grains adopted by Kataoka et al. (2015a), which are spheres with a 
mixture of silicate (8%), water ice (62%) and organics (30%). All 
fractional abundances are in volume and are taken from Pollack et 
al. (1994). We assume a canonical gas-to-dust mass ratio of 100, 
and use the Mie theory to calculate the absorption and scattering 
opacities (Bohren & Huffman 1983). The inferred scattering opac¬ 
ity corresponds to a grain radius a - 31 jjm for grains of a single 
size. For the MRN-type power-law distribution n(a) oc a~^'^ (Mathis 
et al. 1977), we obtain a maximum grain size of amax = 72 fim. The 
increase of this maximum over the single size case comes from av¬ 
eraging over the grain size. In both cases, the dimensionless param¬ 
eter X = Inal A 1, so that the Rayleigh limit used for treating the 
scattering in the previous sections is self-consistent (see Fig. 6 be¬ 
low). The maximum size inferred for the grains responsible for the 
scattered dust emission in the HL Tau disk is much larger than that 
of the grains in the diffuse interstellar medium. This is consistent 
with other lines of evidence for grain growth in disk environments 
(e.g., Perez et al. 2012; Alecian 2013; Testi et al. 2014). 

We note that Kataoka et al. (2015b) independently modeled 
the HL Tau disk polarization using dust scattering through Monte 
Carlo radiative transfer simulations. They obtained disk polariza¬ 
tion patterns that are very similar to ours. They inferred a maximum 
grain size that ranges from 70 jum to 350 jum, which is broadly con¬ 
sistent with our value of 72 fim. 

In summary, to reproduce the ~ 1% polarization fraction ob¬ 
served in the disk of HL Tau through dust scattering, the grains 
must have grown to tens of microns (the exact value depends on 
the assumed grain size distribution and composition). However, this 
picture is complicated by the opacity spectral index p inferred for 
HL Tau, as we discuss next. 

4.2 Opacity Spectral Index p and the Need for Larger 

Grains? 

Kwon et al. (2011) obtained a best-fit value 0.73 for the spectral 
index p of the dust opacity Kabs ^ ^ for tho HL Tau disk based 
mostly on CARMA observations at 1.3 and 2.7 mm. It is in agree¬ 
ment with the spatially averaged value obtained from ALMA ob¬ 
servations from 0.87 to 2.9 mm (ALMA Partnership: Brogan et 
al. 2015). This value is significantly lower than the typical ISM 
value of p ~ 1.5-2. The difference is usually taken as evidence 
for grain growth to millimeter size or larger (Testi et al. 2014), 
although other interpretations are possible. For example, Ricci et 
al. (2012) showed that a value of ~ 1 or lower can be obtained 
without mm/cm sized grains if part of the disk is optically thick. 
Some support for this possibility is provided by the spatially re¬ 
solved distribution of p derived from the ALMA data, which shows 
~ 0 indicative of optically thick emission at the central contin¬ 
uum peak and two rings (B1 and B6, ALMA Partnership: Brogan 
et al. 2015, see their Fig. 3). Another possibility is that the index 
P is sensitive to not only the size but also the shape of the grains. 
Indeed, Verhoeff et al. (2011) was able to reproduce the spectral 
energy distribution (SED) of the disk of HD 142527 (with ~ 1 
in the millimeter regime) with irregular grains of sizes up to only 
2.5 yum; the grain shape was treated with the distribution of hollow 
spheres Min et al. (2005). The grains inferred in our model of dust 
scattering-induced polarization for the HL Tau disk have a signif- 
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Figure 6. Comparison of the distributions of the degree of polarization (see 
text for definition) as a function of scattering angle for three cases: the 
Rayleigh limit (blue line) and two cases with MRN-like power-law grain 
size distribution with amax = 72 jum (green) and 4 mm (red). The Rayleigh 
and amax = 72 jum cases are almost indistinguishable. The negative “peak” 
around 130° in the amax = 4 mm case is an example of the so-called “po¬ 
larization reversal,” which may provide a way to probe large, mm/cm sized, 
grains through scattering-induced polarization. 


icantly larger maximum size (of order tens of microns). They may 
still be able to reproduce the observed (averaged) opacity spectral 
index of j3 ^ 0.73 if the grains are irregular and/or part of the disk 
is optically thick. Detailed exploration of this possibility is beyond 
the scope of the present work. 

If large, mm/cm sized, grains are responsible for the relatively 
low value of p observed in the HL Tau disk, it is natural to ask 
whether they can produce a polarization pattern that matches the 
observed one through scattering. It is unlikely, because the key to 
producing the observed pattern is the polarization degree of the 
scattered light peaking near 90° (as in the Rayleigh limit), and this 
requirement is not satisfied for mm/cm sized grains. For example, 
for the grain model adopted by Kataoka et al. (2015a), the polar¬ 
ization degree (defined as the ratio of the two elements in the scat¬ 
tering matrix, -Z 12 /Z 11 , which is essentially the polarization frac¬ 
tion but can be either positive or negative) is nearly zero at 0.87 
mm for all scattering angles except around 135°, where it reaches 
a (negative) “peak” value of ~ -0.2 for Umax - 1 mm and 1 cm 
(see the right panel of their Fig. 2). The negative value is known 
as the polarization reversal (e.g., Murakawa 2010; Kirchschlager & 
Wolf 2014) which, together with the shift of the polarization “peak” 
away from 90°, is expected to produce a polarization pattern very 
different from the Rayleigh scattering case. 

As an illustration, we repeat the computation of the scattering- 
induced polarization at T = 1.3 mm in § 3, but with an MRN-type 
power-law size distribution up to Umax - 4 mm (instead of 72 jjm), 
using the dust model of Kataoka et al. (2015a) and Mie theory. 
The maximum grain size is chosen such that Umax ~ 3/1, which is 
roughly the minimum value required to yield an opacity spectral 
index of ~ 1 according to Draine (2006). The distribution of the 
polarization degree with scattering angle in this case is shown in 
Fig. 6. It is very similar to that obtained by Kataoka et al. at 0.87 
mm, except that the “peak” is slightly lower (-0.17) and is shifted 
to a slightly smaller angle of ~ 130°. 


1 10 '^ 


Figure 7. Scattering-induced polarization by large grains. As in Fig. 5, plot¬ 
ted are the polarized intensity (color map) and polarization vectors (line 
segments, with length proportional to the polarization fraction). Note the 
strong asymmetric with respect to the major axis in both the polarized in¬ 
tensity and the polarization vectors. The polarization along the major axis 
in the central region is due to polarization reversal, which may be a robust 
indicator of scattering by large, mm/cm-sized, grains. The near side of the 
disk is on the right. 


In Fig. 7, we plot the distribution of the polarized inten¬ 
sity together with polarization vectors for the large grain case of 
(^max - 4 mm. There are several features that are worth noting. 
First, unlike the Rayleigh scattering case, the polarized intensity 
is no longer symmetric with respect to the major axis. This is be¬ 
cause large, mm/cm sized, grains preferentially scatter light in the 
forward direction (e.g., Bohren & Huffman 1983), making the side 
of the disk closer to the observer (the right half) brighter. The po¬ 
larization fraction is, however, higher on the far side (especially 
toward the outer part of the disk) because the polarization degree 
of the scattered light is higher for backward scattering than for for¬ 
ward scattering (see Fig. 6). The most striking difference between 
this case and the Rayleigh scattering case shown in Fig. 5 lies in 
the polarization direction. The difference comes from the polar¬ 
ization reversal in the large grain case, which yields an intrinsic 
(or face-on) polarization direction in the radial (as opposed to az¬ 
imuthal) direction and an inclination-induced polarization along the 
major (rather than minor) axis. The interplay between the intrin¬ 
sic and inclination-induced polarization leads to polarization direc¬ 
tions in the region of high polarized intensity (the most easily ob¬ 
servable part) completely different from those observed in HL Tau 
(see Fig. 1). 

We are therefore left with an interesting conundrum. The po¬ 
larization pattern in the HL Tau disk is suggestive of Rayleigh scat¬ 
tering by relatively small dust grains (although still much larger 
than the typical ISM grains), but such grains may have difficulty 
reproducing the observed opacity spectral index(< 1). The index 
can be reproduced more easily with larger, mm/cm sized, grains, 
but it is difficult to generate the observed polarization pattern with 
such grains through scattering. It is conceivable that there are two 
populations of dust grains, with one responsible for polarization, 
the other for p. The two populations do not have to be located 
CO- spatially in the disk; for example, large grains responsible for 
the bulk of the unpolarized continuum (and thus p) may have set- 
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tied close to the midplane, whereas smaller grains that dominate 
the polarized millimeter radiation may remain floating higher up 
above the midplane (e.g., Dullemond & Dominik 2004; Tanaka et 
al. 2005; Balsara et al. 2009 ). If this speculation turns out correct, 
polarized emission in millimeter would provide a powerful probe of 
not only grain growth, but also the expected vertical stratiflcation 
of grain sizes, especially in conjunction with observations of opti- 
cal/TR polarization, which probe even smaller, micron-sized, grains 
that are higher up still above the disk midplane. 


5 OTHER SOURCES AND FUTURE DIRECTIONS 

As mentioned in § 1, spatially resolved polarized sub-mm/mm 
emission was observed in two other sources besides HL Tau: IRAS 
16293-2422B (Rao et al. 2014) and L1527 (Segura-Cox et al. 
2015), using SMA and CARMA respectively. We discuss whether 
the polarization in these two sources is compatible with an origin 
in dust scattering. 

The LI527 disk is nearly edge-on, with an inclination angle 
i ^ 90°. Its 1.3 mm dust emission was observed to be polarized 
at ~ 2.5% level, with a direction roughly perpendicular to the disk 
(i.e., along the minor axis). This pattern is consistent with that ex¬ 
pected from dust grains aligned by a predominantly toroidal mag¬ 
netic held (Segura-Cox et al. 2015). It is also the pattern expected 
of scattering by relatively small dust grains in the Rayleigh limit 
because the polarization induced by disk inclination is along the 
minor axis (as illustrated vividly in the lower right panel of Fig. 2). 
Indeed, the fraction of the scattered radiation that is polarized due 
to disk inclination increases with the inclination angle, reaching 
a maximum value of 1/3 for edge-on disks, as we showed ana¬ 
lytically in § 2.2 (see eq. (18)). This trend makes it more likely 
for L1527-like disks to show scattering-induced polarization than 
face-on disks. Indeed, the median polarization fraction of L1527 
is the highest among the three sources with spatially resolved sub- 
mm/mm polarized emission so far. However, whether such polar¬ 
ization can actually be observed or not depends on the polarized 
intensity, which in turn depends on the total intensity and the polar¬ 
ization fraction (relative to the total intensity). The latter is sensitive 
to the scattering opacity Ksca, which depends on the dust properties, 
especially the grain size, which can vary greatly from one source to 
another. To produce the observed polarization fraction of ~ 2.5% 
in L1527 through scattering, the grains must have grown well be¬ 
yond fim size; otherwise, the scattering opacity would be too small. 
The scattering cannot be dominated by large, mm/cm sized grains 
either. Such grains would produce polarization along, rather than 
perpendicular to, the edge-on disk because of polarization reversal 
(see § 4.2 and Figs. 6 and 7). Large cm-sized grains are inferred in 
the LI527 disk from the small opacity spectral index ~ 0 (Tobin 
et al. 2013). If the observed polarization is due to scattering by sub- 
mm sized grains, then two grain populations may again be needed, 
as in the HL Tau case discussed in § 4.2. 

The disk in IRAS 16293-2422B appears to be nearly face-on 
(Rodriguez et al. 2005; Zapata et al. 2013). Its 0.88 mm emission 
was observed to be polarized at ~ 1.5% level, with the polarization 
directions showing a swirling pattern that is neither strictly radial 
nor purely azimuthal (Rao et al. 2014). Rao et al. showed that the 
pattern is broadly consistent with the polarized emission expected 
from grains aligned by a magnetic held that is warped into a spiral 
conflguration by disk rotation, although the rotation is hard to mea¬ 
sure directly because of the face-on orientation. It is inconsistent 
with the simplest dust scattering model for axisymmetric face-on 


disks, where the polarization directions are expected to be perfectly 
azimuthal, as illustrated in the upper left panel of Fig. 2. However, 
the observed polarized intensity is arc-shaped which, for face-on 
systems, requires intrinsically non-axisymmetric disk models, such 
as those constructed by Kataoka et al. (2015a). Indeed, the observed 
intensity distribution is reminiscent of the lopsided ring models of 
Kataoka et al. (see their Fig. 6 and 7), although the polarization di¬ 
rection in the models changes sharply from radial in the inner part 
of the ring to azimuthal in the outer part. It is conceivable that dust 
scattering models with more complicated disk structures, such as 
spiral arms, may match the observation better, but this remains to 
be demonstrated. One worry is that the dust emission in this source 
may be optically thick at sub-mm wavelengths (Zapata et al. 2013), 
which would reduce the degree of anisotropy in the unpolarized 
radiation to be scattered, and thus the degree of polarization in the 
scattered radiation (Kataoka et al. 2015a). Another is that the polar¬ 
ized emission detected in this deeply embedded source may be con¬ 
taminated by the protostellar envelope. Higher resolution ALMA 
polarization observations should become available in the near fu¬ 
ture. They will provide more stringent tests of the dust scattering 
model of polarized mm/sub-mm emission. 

As stressed by Kataoka et al. (2015a), polarized radiation 
at mm/sub-mm wavelengths provides a powerful probe of grain 
growth, if it is produced by dust scattering. A robust prediction 
of the scattering model is that large, mm/cm-sized, grains should 
produce millimeter polarization along the major axis of an inclined 
disk due to polarization reversal, especially in the central region 
where the intrinsic (face-on) polarization is expected to be weak 
and the observed polarization pattern is more easily dominated by 
that induced by disk inclination (see Fig. 7 for an illustration). This 
effect should be searched for with high resolution ALMA observa¬ 
tions, especially in high-inclination systems. Another way to probe 
large grains is to observe polarization at longer, centimeter, wave¬ 
lengths using, for example, JVLA and the future SKA. At such 
wavelengths, the scattering by mm-sized grains would still be in the 
Rayleigh regime, with a high degree of polarization of the scattered 
light. Given the sensitive dependence of the scattering opacity Ksca 
(which controls the polarization fraction of the total intensity) on 
the grain size relative to the wavelength, it is important to carry out 
high-resolution polarization observations at multiple wavelengths 
to determine whether the polarization is indeed due to dust scat¬ 
tering and, if yes, to constrain the grain size distribution. With the 
polarization capability of JVLA, and that of ALMA coming online 
soon, the prospect for using sub-mm/mm/cm polarization to probe 
grain properties in disks is bright on the observational side. 

On the theoretical side, much work remains to be done. In this 
paper, we have limited our treatment to the simplest case of opti¬ 
cally and geometrically thin disk, to bring out the essential features 
of the disk inclination-induced polarization transparently through 
a semi-analytic model. As noted earlier, there is indication from 
the spatial distribution of the opacity spectral index (5 that part of 
the HL Tau disk is optically thick at sub-mm/mm wavelength, in¬ 
cluding the central continuum peak, B1 and B6 rings, at radii of 
~ 0, 20, and 81 AU, respectively (ALMA Partnership: Brogan et 
al. 2015). Since the B6 ring is rather narrow, and the bulk of polar¬ 
ized emission is detected interior to it, it should not affect our HL 
Tau model much. The continuum peak and B1 ring are not resolved 
by the polarization observation, but they can potentially lower the 
polarization fraction in the central pixel relative to the outer part, 
which is expected to bring the dust scattering model presented in 
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§ 3 and Fig. 5 into closer agreement with observation'*. The effects 
of optically thick regions, as well as substructures on the disk such 
as rings and gaps, need to be quantified in future calculations in 
order to compare with higher resolution ALMA polarization obser¬ 
vations that should become available soon. 

Another future improvement is to relax the thin disk approx¬ 
imation. While the approximation is adequate for large grains that 
have settled close to the disk midplane, it would be less so for 
smaller grains that remain higher up above the midplane. It is likely 
that there is a gradient in dust grain concentration and size distribu¬ 
tion in all three (radial, vertical and azimuthal) directions through 
grain growth, inward drift, vertical settling, and trapping (see Testi 
et al. 2014 for a recent review, and Perez et al. 2012 for an example 
of gradient in grain size). Such gradients should be taken into ac¬ 
count in more complete models of scattering-induced polarization, 
perhaps using 3D radiative transfer codes such as RADMC3D. As 
noted earlier, in HL Tau and LI527, the observed opacity spectral 
index p and polarization pattern in millimeter, if originating from 
dust scattering, appear to require two grain populations of different 
size distributions. Whether they can arise naturally through grain 
evolution in the disk warrants further investigation. In addition, a 
complete model of disk polarization at sub-mm/mm/cm will need 
to include both the polarized emission by aligned, non-spherical 
grains (particularly grains of tens of microns in size or larger, al¬ 
though it is unclear whether such grains can be aligned with respect 
to the magnetic field or not, as discussed in § 1) and dust scattering, 
and possible interplay between the two. 


6 CONCLUSIONS 

Motivated by the recent spatially resolved millimeter polarization 
observations of the HL Tau disk, we have developed a simple semi- 
analytic model for the dust scattering-induced polarization in the 
limit of optically and geometrically thin disk and Rayleigh scatter¬ 
ing, with an emphasis on the effects of the disk inclination to the 
line of sight. The main results are summarized as follows: 

1. We developed an efficient approximate method for comput¬ 
ing disk polarization from dust scattering by dividing the source 
region of the millimeter radiation to be scattered at a location in¬ 
side the disk into two conceptually distinct parts: a near-field re¬ 
gion centered on the location with a size comparable to the local 
dust scale-height, and a far-field region outside. Radiation from the 
near-field region is more or less isotropic, and does not contribute 
significantly to the polarization of the scattered light. Radiation 
from the far-field region is concentrated toward the disk plane. It 
is strongly polarized after scattering in an inclined disk. The polar¬ 
ization fraction of the scattered light increases with the inclination 
angle, reaching a maximum value of 1/3 for edge-on disks if the 
incoming radiation to be scattered is azimuthally isotropic in the 
disk plane (eq. (18)). The polarization induced by disk inclination 
is parallel to the minor axis. It can easily dominate the intrinsic 
polarization of the disk in the face-on view (see Fig. 2). 

2. We developed a simple model for the polarization of the HL 
Tau disk, based on the(Kwon et al. 2011) model of disk physical 
structure and polarization induced by a disk inclination of 45° (see 
Fig. 5). The model naturally reproduces two main features of HL 
Tau: (1) the region of high polarized intensity is elongated along 

The polarization fraction of millimeter emission in the central region can 
also be lowered if the grains there grow to centimeter sizes, which can con¬ 
tribute significantly to the total but little to the polarized intensity. 


the major axis, and (2) the polarization vectors in this region are 
roughly parallel to the minor axis. Both are the consequences of a 
simple geometric effect: only the radiation propagating along the 
major axis of a tilted disk would be scattered by 90° to reach the 
observer and be maximally polarized, with a polarization direction 
along the minor axis in the plane of the sky. The broad agreement 
is robust because it does not depend on the detailed properties of 
dust grains (which are uncertain) as long as the scattering is in the 
Rayleigh limit. It provides support for the millimeter polarization 
observed in this particular case originating at least in part from dust 
scattering, although polarized emission from magnetically aligned 
dust grains cannot be ruled out, especially if the disk field is more 
complex than toroidal. 

3. For the other two cases with observed mm/sub-mm polar¬ 
ization, L1527 and IRAS 16293-2242B, the situation is less clear. 
The observed polarization pattern in the nearly edge-on disk of 
LI527 is compatible with that expected of either toroidal field- 
aligned grains or dust scattering. The pattern observed in the pos¬ 
sibly face-on disk of IRAS 16293-2242B is more consistent with 
that expected of grains aligned by a rotationally warped mag¬ 
netic field than with the simplest case of dust scattering. Higher 
resolution observations of more disks with different inclination 
angles are needed to better differentiate the grain-alignment and 
dust-scattering models. The observational situation should improve 
drastically with ALMA and JVLA. 

4. To reproduce the polarization fraction of ~ 1% observed at 
1.3 mm in the HL Tau disk, a maximum size of tens of microns is 
needed for the scattering grains. Such grains are generally thought 
to be too small to produce an opacity spectral index p of order 1 or 
less that is observed in HL Tau and other sources; larger, mm/cm 
sized grains may be needed. However, mm/cm sized grains tend to 
produce polarization parallel (rather than orthogonal) to the major 
axis due to polarization reversal (see Figs. 6 and 7), which is not 
observed in HL Tau; nevertheless, this pattern should be searched 
for in other sources as a robust indicator for large grains. In any 
case, the dust scattering model for polarization and the relatively 
low p produce an interesting conundrum that needs to be resolved 
in the future, perhaps with more complete models that include grain 
evolution and 3D radiative transfer, as well as polarized direct emis¬ 
sion from aligned grains. Such models, together with the high res¬ 
olution ALMA/JVLA polarization observations that will soon be¬ 
come available, should make it possible to disentangle the contribu¬ 
tions of grain alignment and dust scattering to the disk polarization, 
which is needed in order to provide robust constraints on the mag¬ 
netic field that is generally thought to be crucial to the dynamics 
and evolution of protoplanetary disks and/or the grain growth that 
may eventually lead to planet formation. 
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